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Summary 

A method of deriving nondimens ional equations and indentifying the 
fundamental parameters associated with bifurcation buckling of shallow 
shells subjected to combined loads is presented. More specifically, 
analysis is presented for symmetrically laminated doubly-curved shells that 
exhibit both membrane and bending anisotropy. First, equations for 
nonlinear deformations of thin elastic shallow shells are presented, and 
buckling equations are derived following the method of adjacent equilibrium 
states. Next, the procedure and rationale used to obtain useful 
nondimensional forms of the transverse equilibrium and compatibility 
equations for buckling are presented. Fundamental parameters are identified 
that represent the importance of both membrane and bending orthotropy and 
anisotropy on the results. Moreover, generalizations of the well-known 
Batdorf Z parameter for symmetrically laminated shells with full anisotropy 
are presented. Using the nondimensional analysis, generalized forms of 
Donnell's and Batdorf' s equations for shell buckling are also presented, and 
the shell boundary conditions and approximate solution methods of the 
boundary -value problem are briefly discussed. 

Results obtained from a Bubnov-Galerkin solution of a representative 
example problem are also presented. The results demonstrate the advantages 
of formulating the analysis in terms of nondimensional parameters and using 
these parameters for parametric studies. The results specifically show that 
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shells with positive Gaussian curvature are much more shear buckling 
resistant than corresponding flat plates and shells with negative and zero 
Gaussian curvature. In addition, the results indicate that the importance 
of bending anisotropy on shear buckling resistance Is influenced by shell 
curvature . 


Introduction 

Buckling behavior of laminated composite plates and shells has received 
renewed interest in recent years due to the search for ways to exploit 
anisotropy in the design of aerospace vehicles. The present study focuses 
on buckling analysis of symmetrically laminated doubly-curved shallow 
shells. These shells are candidates for aircraft applications such as wing 
cover panels, empennage and fuselage skins, engine cowlings, and wing-to- 
body fairings. In addition, these shells have potential spacecraft 
applications such as liquid fuel tankage, pressure bulkheads, missle nose 
cones, and payload modules. Understanding the fundamental parameters that 
affect the performance of these shell structures, as well as the importance 
of anisotropy, is a key ingredient of their design, and ultimately their 
use. Moreover, these parameters can provide valuable Information useful in 
developing scaling laws for structural testing of plates and shells. 

Symmetrically laminated shells exhibit anisotropy in the form of 
material- induced coupling between pure bending and twisting deformations 
that is commonly exhibited by symmetrically laminated plates undergoing 
bending deformations. In addition, they exhibit anisotropy in the form of 
coupling between pure biaxial stretching and membrane shearing. The 
anisotropy is manifested in shell buckling theory by the presence of odd- 
combination mixed partial derivatives in the partial differential equations 
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governing buckling and the natural or force boundary conditions. The 
presence of the these derivative terms prevent simple closed form solutions 
from being obtained. Studies that assess the importance of anisotropy on 
buckling behavior of symmetrically laminated plates are presented in 
references 1 and 2. In reference 1, nondimensionalization of the partial 
differential equation governing plate buckling is presented and 
nondimens ional parameters are identified that provide physical insight into 
plate buckling behavior. The results presented in reference 1 suggest the 
potential for using nondimens ional parameters to characterize plate buckling 
behavior, assess the importance of anisotropy, and aid in their preliminary 
design. The results presented in the present paper build upon the results 
presented in reference 1. 

A major objective of the present paper is to present nondimens ional 
parameters that will aid in, and simplify, buckling analysis and preliminary 
design of laminated composite shallow shells. It is a well-known fact that 
shells with a high degree of curvature are sensitive to small imperfections 
in their geometry under certain loading conditions, and that this 
imperfection sensitivity leads to collapse loads often substantially below a 
predicted bifurcation buckling load [3], However, there is a class of shell 
problems for which imperfection sensitivity is minimal under certain loading 
conditions and the addition of a slight amount of curvature has a positive 
effect on increasing the buckling resistance of a flat plate [4] . In this 
case, results obtained from a bifurcation buckling analysis can be used to 
obtain credible estimates of the collapse load. It is for this class of 
shallow shells that the nondimens ional parameters and equations presented in 
the present paper have been derived. A significant example of a 
nondimens ional parameter that has seen wide use in the buckling analysis and 
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design of isotropic cylindrical shells is the well-known Batdorf Z or 
curvature parameter [5,6]. 

Nondimens ional parameters can play a key role in the preliminary design 
of aerospace vehicles. An important consideration is that, in the initial 
preliminary design phase, a structural designer is usually under severe time 
constraints and often prefers to have information available in a handy chart 
form. Nondimensional parameters permit fundamental results to be presented 
as a series of curves, on one or more plots, that cover the complete range 
of shell dimensions, loading combinations, boundary conditions, and material 
properties. In addition, the curves also furnish the designer with an 
overall indication of the sensitivity of the structural response to changes 
in geometry, loading conditions, boundary conditions, or material 
properties. Often in preparing design charts of this nature, a special 
purpose analysis is preferred over a general purpose analysis due to the 
cost and effort involved in generating results encompassing a wide range of 
design parameters. Examples of design charts for buckling and postbuckling 
of orthotropic plates that use nondimensional parameters are presented in 
references 7 and 8. Another very useful benefit of nondimensional 
parameters is that they can be used to greatly simplify the equations and 
results when performing analysis and to provide insight into the order and 
importance of various terms appearing in the equations. 

Two major objectives of the present paper are to present a method of 
deriving nondimensional equations for doubly- curved shallow shells subjected 
to combined loads, and to indentify the fundamental parameters associated 
with bifurcation buckling of these shells. The paper begins with a 
presentation of the equations for nonlinear deformations of symmetrically 
laminated elastic shallow thin shells. Next, the equations for buckling are 
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presented using the method of adjacent equilibrium states, and are cast into 
a form suitable for nondimensionalization. The paper then indicates the 
steps required, and the rationale used, to obtain nondimens ional forms of 
the transverse equilibrium and compatibility equations for buckling. 
Fundamental parameters are then identified that can be used to represent the 
importance of membrane and bending orthotropy and anisotropy on shell 
behavior. Parameters are also presented that are generalizations of the 
Batdorf Z parameter for symmetrically laminated shells with full anisotropy. 

After obtaining the nondimens ional buckling equations and identifying 
the nondimens ional parameters, generalized forms of Donnell's and Batdorf' s 
equations for shell buckling are presented. These generalized equations 
contain the effects of both laminated material properties of the shell as 
well as double curvature. Finally, a brief discussion of shell boundary 
conditions and approximate solution methods for the shell buckling problem 
is presented. 

The analysis presented in the present paper was inspired by the work 
presented in reference 8. For this reason, and for many useful discussions, 
the author would like to dedicate this paper to the late Dr. Manuel Stein of 
NASA Langley Research Center. 

Equations for Nonlinear Deformations 

The basic equations for doubly-curved shallow shells of general shape 
and uniform thickness t are presented in this section in terms of the 
orthogonal lines-of -curvature curvilinear coordinates shown in 

figure 1. The equations presented consist of the strain-displacement 
relations, nonlinear equilibrium equations, compatibility equations, and 
constitutive equations for symmetrically laminated shells. 
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Strain-Displacement Relations 

The nonlinear strain-displacement relations used herein to describe the 
deformation of a shallow shell of general shape are the relations of 
Donnell -Mushtari -Vlasov classical thin-shell theory. These relations are 
well known (for example see reference 6, pp. 190-197), and are given by 
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The symbols and u 2 are the displacements of the shell reference surface 
(in the tangent plane) in the £ and $ 2 directions, respectively, and the 


symbol w is the displacement along the direction normal to the shell 
surface. The coordinate f corresponds to a distance along an axis normal to 
the shell surface, at a given point, as shown in figure 1. The terms and 

R 2 represent the principal radii of curvature of the reference surface along 

the and £ 2 coordinate directions, respectively, and the symbols and A 2 

are the coefficients of the first fundamental form of the surface given by 

(dS) 2 - (A 1 d^ 1 ) 2 + (A 2 d£ 2 ) 2 (10) 


where dS is the differential arc-length between to neighboring points on the 
surface. The terms £°, £°, and y° 2 are the membrane strains, and the terms 

K °, «° and k ° u , when multiplied by f are the strains associated with 

o o 

changes in curvature of the shell reference surface. The terms and y 23 

are the transverse shearing strains and are assumed to be zero -valued in 
accordance with the assumptions of classical first-approximation thin-shell 
theory. Enforcing this condition yields expressions for the surface 
rotations in term of the normal deflection w, i.e. , 
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Nonlinear Equilibrium Equations 

The nonlinear equilibrium equations of Donnell -Mushtari -Vlasov shallow 
shell theory are obtained by enforcing equilibrium of a differential shell 
element in its deformed configuration. The equations resulting from this 
process are given as follows. The equation corresponding to the summation 
of forces in the direction is given by 
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where „ N 2 , and N^ 2 are the membrane stress resultants and q^ is an 

applied membrane surface traction. Similarly, the equation corresponding to 
the summation of forces in the £ 2 direction is given by 
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where q^ is also an applied membrane surface traction. The equation 
corresponding to the summation of forces in the f direction is given by 
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where P denotes the contribution of the membrane forces to the transverse 
m 

(normal to the surface) equilibrium that is given by 
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(15b) 


+ f? 2 [N 12 A l ( A 1 f? 1 ) + N 2 A l ( A 2 f? 2 )] 


The symbols and Q 2 denote the transverse shearing stress resultants. The 
equation corresponding to the summation of moments about the £^-axis is 
given by 
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where M^, Mg, and M^g are the bending stress resultants. Similarly, 
equation corresponding to the summation of moments about the fg-axis 
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A more convenient form of the nonlinear equilibrium equations, for the 
purpose of defining nondimensional parameters, is obtained by requiring the 
curvilinear coordinates and £g to be surface coordinates with units of 


length. A new set of curvilinear coordinates s^ and Sg are introduced such 


that 

(dS) - (ds 1 ) + (ds 2 ) (18) 

For this set of coordinates, A^ — 1 and A 2 = 1. The nonlinear equilibrium 


equations simplify to 
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A similar set of equations derived for postbuckling analysis of isotropic 
shallow shells subjected to thermal loads is presented in reference 9. 

The equilibrium equations given in equation (19) are reduced to a set 
of three equilibrium equations as follows. First, the expression for 


given by equation (19d) is rearranged, by differentiating, to give 
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Substituting equations (19a) and (19b) into equation (20) gives 
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Next, eliminating the transverse shear stress resultants from equation (19c) 
using equations (19e) and (19f) gives 
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Equations (19a), (19b), and (22) are the nonlinear equilibrium equations 
used in this paper to obtain the equations for buckling of shallow shells. 


Compatibility Equation 

The compatibility equation for the nonlinear boundary- value problem is 
obtained by eliminating u^ and u^ from equations (2), (3), and (4). It is 


convenient, 
through (4) 


for the purpose of the present study, to express equations (2) 
in term of the length coordinates; i.e., 
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For shallow shells, for which the curvatures are mildly varying, the 
following approximation is used to simplify the compatibility equation; 
i.e. , 


2 

/W_ 

2 '■R 

a Sl 2 


) 


oc 


1 <L_w 

R 2 

R 2 as x 


(27a) 


11 


(27b) 


2 

*-?<£-> 

Ss 2 1 


1 3 V 

D 2 

1 5s 2 


Using this approximation, the compatibility equation becomes 
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Constitutive Equations 

The shells considered in this paper are symmetrically laminated and 
exhibit both anisotropic membrane and bending behavior. The corresponding 

constitutive equations are given by 
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where A.., A, „ , A.., and A, , are the orthotropic membrane stiffnesses; 

11 12 22 oo 

Aj^ and A 2 g are the anisotropic membrane stiffnesses; D.^, °22’ anc * °66 

are the orthotropic bending stiffnesses; and and are the anisotropic 
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bending stiffnesses of classical laminated thin-shell theory (for example 
see reference 10) . 


Buckling Equations 

The equations governing buckling are derived using the method of 
adjacent equilibrium described in reference 6 (see pp. 201-202). Prior to 
buckling, the shell is assumed to be in a primary equilibrium state given by 

the displacement field (u*\ u^, w P ) . Near the point of buckling, there is 

assumed to exist an adjacent equilibrium state (Su^ F 6 w) sufficiently 

close to the primary equilibrium state such that Su^ f ^u^, and 6 w are 

infinitesimal quantities. Moreover, prebuckling displacements normal to the 

shell reference surface are presumed to be negligible; i.e., * 0. Using 

this method, the linearized bifurcation buckling equations are obtained by 
substituting 

u^ -* u^ + 8u ^ (31a) 

u 2 - U P + 6u 2 (31b) 

w -+ 0 + 6w (31c) 

into the basic equations of the nonlinear boundary -value problem given 
previously. Associated with the substitutions defined by equations (31) are 
variations in the strains and changes in curvature. In addition, increments 
in strains and changes in curvature produce increments in the stress 
resultants via the constitutive equations. The relationships for the 
strains and changes in curvature are given by 
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and where the superscript p denotes the primary equilibrium state, 
relationships for the stress resultants are given by 
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and result from substituting equations (32) into the constitutive equations 
(29) and (30). The superscript p again denotes the primary equilibrium 
state. In equations (32) and (33) terms of quadratic degree and higher in 
5u^, 5^, and 5 w have been omitted in accordance with the assumption that 

the variations are infinitesimal in size. The specific form of the 
constitutive equations are given by 
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The form of the buckling equations used herein are obtained by 
substituting equations (31) through (34) into equations (19a), (19b), (22), 
and (28), Neglecting terms greater than linear in the displacement, strain, 
changes of curvature, and stress resultant increments yields the following 
equations of the adjacent equilibrium state (buckling equations) 
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where the superscript p on the prebuckling stress resultants has been 
dropped for simplicity. An additional operation used, but not explicitly 
shown, to obtain the buckling equations given above is the enforcement of 
the requirement that the primary equilibrium state satisfies the equilibrium 
equations and compatibility equation of the nonlinear boundary -value 
problem. 
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Nondimensionalization of the Buckling Equations and Parameters 
One objective of this paper is to determine the nondimens ional 
parameters governing the buckling behavior of shallow thin elastic shells. 
This objective is accomplished following the procedure presented by Batdorf 
in references 11 through 13 for isotropic curved plates, by Stein in 
reference 8 for flat laminated orthotropic plates, and by the author in 
reference 1 for flat symmetrically laminated plates. The rationale used in 
this paper for deriving the nondimens ional equations and selecting the 
nondimens ional parameters follows, to a great extent, the work presented by 
Stein [8] . This rationale is to make the the field variables and their 
derivatives of order one, minimize the number of independent parameters 
required to characterize the behavior, and to avoid introducing preferential 
direction into the nondimensional equations. 

The first step in the nondimensionalization procedure is to formulate 
the boundary-value problem in its simplest form which corresponds to two 
coupled homogeneous linear partial differential equations. This task is 

accomplished by introducing a stress resultant function defined by 
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This stress resultant function (or simply stress function) satisfies 
equations (36a) and (36b) identically. Upon eliminating these two equations 
from the boundary -value problem, the equation guaranteeing compatibility of 
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the buckling strains must be satisfied* This buckling compatibility 
equation is given by equation (36d). To get the buckling compatibility 
equation into a solvable form, the inverted form of the constitutive 
equations are used to express the buckling strains in terms of the stress 
function; i . e . , 
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where equations (37) have been substituted into the usual constitutive 
equations involving the membrane stress resultants, and where a^ , a l2* a 22* 


a. . , and a, r denote the coefficients of the inverted form of the 
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membrane constitutive equations. Substituting equations (38) into the 
buckling compatibility equation (36d) and simplifying gives 
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Following the analysis presented in reference 8, the following 
nondimensional coordinates are used 
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where and I , are characteristic dimensions of the shell shown in figure 
1. Substituting these nondimens ional coordinates into equation (39) and 
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yields 


2 LJ± ^ o 3-1$- A o d 6$ . 9 _jn A , 1 

a 7" + 2a 7 3 +2/i 2 2 + 2 3 + 2 4 

dz- dz^dz n dz,dz 0 m dz,dz„ a dz„ 

1 12 12 12m 2 


where 


2 

L, 2 


2 

L, 2 


7 a ii a 22 1 r 2 az * 


"2 d 6w "1 3_£w, 

° R 2 J 


1 3z, 


m 


12, , .1/4 

L 1 (a 22 /a ll ) 


2a 12 + a 66 


2/ a U a 22 


26 


T nT * 3 .1/4 

( a H a 22 ) 


and 


16 


m 


, 3 a/4 

(a ll a 22 ) 


(41) 


(42a) 


(42b) 


(42c) 


(42d) 


The next step in the nondimensionalization involves the buckling 
equation given by equation (36c). To express this equation in terms of the 
buckling displacement 5w and stress function 5$, the equations defining the 
stress function and the constitutive equations given by equations (35d) 
through (35f) are used. Substituting equations (37) and (35d) through (35f) 
into equation (36c) and simplifying gives 
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4 

d_6w 


4 

3 3w 


D 11 < + 4D 16 3 
os^ as^os^ 


+ 2(D 12 + 2D 66 )_ ^T + - 2 - 45£ J + D 

as ds 2 ds l ds 2 ds 2 


2 2 2 2 

i a 5$ i a 5$ aaw 3Sw a 5 w a aw 

R ’ + R, > + "l 3s, + q 2 3s - B 1 7^ ' N 2 . » - 2N 

1 3s ^ 2 3s^ 1 2 3s^ 3s2 


2 

3 8 w 

12 


- 0 


(43) 

Substituting the nondimens ional coordinates defined by equation (40) into 


2 t 2 


equation (43) , multiplying through by L^L2/7 d^D 22 yields 


3 S w 


-b . “ + Vb * 


3 5w 


+ 2$ 


dz. 


3z^3z2 


J LJj d- 
2 2 
3z^3z2 


+ 4‘ 


b 3 <5w 1 3 3w 


2 4 


b 3 z-3z 0 a, 3z* 

1 z b z 


3^<$> 


2 3 5$ ~ 3£w 

2 t — 2 + q - + q 

V D 11 D 22 5z 2 R 2^ D 11 D 22 dz l 1 


36w 
2 


2 3 5w 
3z^ 


2 <3 6w 
v o' r 2 

a z 2 


K 10 7r o 2 c 

12 3 6w 


ot. 3 z- 3z- 
b 1 z 


(44) 


where the nondimensional parameters appearing in this equation are given by 

L„ 


»b - ^ ([ V D 22> 1/4 


(45a) 


P T2 + 2D 66 

7 D 11 D 22 


(45b) 


16 


'b 3 1/4 

< D 11 D 22> 7 


(45c) 


D 


26 


b (D 11 D 22 >1/4 


(45d) 
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(45e) 


K,- 


N 1 L 2 


1 * Al D 22 


K„ 


M 2 L 1 


K 


* ^ D 11 D 22 


N 12 L 2 


12 2 3 1/4 

* < d u d 22 ) 7 


(45f) 


(45g) 


q l L l L 2 

y D u D 22 


(45h) 


q 2“ 


q 2 L 2 L l 

y D U D 22 


(451) 


To obtain equations of order one, a new stress resultant function 
defined by 


ST - «*/yD n D 22 


(46) 


is introduced into equations (44) and (41) . Equation (44) becomes 


2 3 £w 

*b a 4 
dz i 


+ 4a b 7 b 3 


L in 


. 0 n 2. gw Jb g gw, 1 

+ 22 +4 ^ 3 + 2 4 


3z^3z2 


2 2 


b dz^dz^ 


“b dz 2 


u h LJ1 . h LM . - Ms - Ms 


R i a 2 

1 gz 2 


R 2 

R 2 3z x 


1 gz. 


*2 gz. 


K 2 OH 2 LM . -^3- . o 

K 1 2 2 2 a. gz,gz 0 

gz L gz 2 b 1 2 


(47) 


Similarly, substituting equation (46) into equation (41) and simplifying 
gives 
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2 


m 


f . „ 

4 + 2q _t 


dz „ 


m m 


a%F 

3 

dz^dz^ 


t 0 a%F 
+ 2 2 


5 4 4 

, 9 _m a <SF 1 a 6F 

+ z a 3+2 4 

az,az 0 m az,az 0 a az„ 
12 1 2 m 2 


2 2 
L 2 I 2 

2 a aw i a gw 

r 2 3z ; + R 1 a z ; 


^ a ll a 22 D ll D 22 


The term given by 


12ya U a 22 D ll D 22j 


1/2 


( 48 ) 


has dimension 1/t, where t is the shell wall thickness. To get the equation 
to a form of order one, a nondimensional displacement 5W is introduced into 
equation (48) ; i . e . , 


SW 


S w 


12'/ a 11 a 22 D 11 D 22 


1/2 


(49) 


The nondimensional displacement <SW defined in equation (49) has character 
similar to 5w/t. Using equation (49), equation (48) simplifies to 


2 £_£F . 9 

“m 4 + 2o _T 


4 4 6 

as? " r ~ 


m 


dz. 


m m 


az 1 az 2 


. . ollD .JUI 1 £_£F 

+ 3 * 2 4 


az 1 az 2 


m az.az, a az 0 
12 m2 


where 


/To i -7 5_£W 


5_5W, 


/12{Z + Z ““} 

3z^ 3 z 2 


1 R 1 7l2 < 


a l 1 a 22 D ii D 2 2 ] 


R 2 7l2 ja 11 a 22 D ii D 2 2 ] 

Equation (47) also simplifies to 


1/4 


(50) 


(51a) 


(51b) 
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a sw 


dz. 


4 + ^ 


a gw 

b' 3 
dz l dz 2 




az x az 2 


b az 1 az 2 


% 3Z 2 


+ /^rz ^ + Z ^1 +5^ + 5^ 

2 ^2 2 J q l dz, q 2 3z 0 

oz. oZ- A- ^ 


„ 2 a%w v 2 a 2 gw , K i2* a 2 *w 

1 3zj 2 tz\ % dz l S *2 


Equation (50) can be expressed in operator form as 


D (5F) - Jl2 D (6 W) - 0 
m c 


where 


( 52 ) 


(53) 


D (6F) - a 2 + 2a 7 ^ £1 * z 

m m dz x m m az x az 2 az x az 2 ra Sz x dz 2 a m dz 2 


2 , * 2“ + \ M 1 (54a) 


and 


MW . Z, ♦ *! H“ 

3z^ 3 z 2 


(54b) 


The operators and D c are referred to herein as the membrane stiffness 

operator rnd the curvature operator, respectively. 

To get equation (52) into the desired form, it is necessary to express 

the nondiriensional stress resultants K^, K 2> K^ 2 , q^, and q 2 (see equations 

(45)) in terms of a loading parameter p. These relationships are given by 

K, 


- -<? 1 g 1 (z 1) z 2 )p 

(55a) 

- -^ 2 g 2 (z 1 ,z 2 )p 

(55b) 

, - Ag (z^z^p 

(55c) 
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(55d) 


^1 “ ^4g 4 ( z i- z 2>P 

q 2 “ / 5 g 5 ( z 1 . z 2 ) P ( 55e ) 


where the minus signs are used to make compression loads produce positive 
values of the critical loading parameter. The parameters / ^ / 3 * , 

and^<- are load factors that indicate the ratio of the nondimens ional 


membrane stress resultants and surface tractions prior to buckling. The 
functions g^(z^ f z^) through g^(z^ t z^) indicate the spatial variation of the 


nondimens ional membrane stress resultants and surface tractions. Using 
these relations, equation (52) is expressed in operator form as 


D b (<5W) + 712 D c (5F) - pK («W) 


(56) 


where the operators are defined by 


D b (fiW) 


4 4 

2 a SVJ . 8 8 w 

-b rr + 4< v b : 

dz ^ dz^dz^ 


. _OSL 

■F ^P 2 2 


fb sLM l 

+ 3 + 2 4 

dz^dz r. b 3z-3z- a, dz 0 

12 1 2 b 2 


(57) 


K (5W) 
g 


* 2 , , 2 /jw 

, Z 2 ) 7T 2 ** ^2®2^ Z 1 ,z 2 ^ 2 

dz^ dz^ 


2 ^3 g 3^ 


z r z 2 ) 


a%w , , ,MH „ .MW 

a Z]L az 2 ' '^4 g 4 (z i ,z 2 ) az 1 '^5 g 5 (z i’ 2 ) az 2 


(58) 


The operators and are referred to herein as the bending stiffness and 

geometric stiffness operators, respectively. 

Equations (53) and (56) constitute the eigenvalue boundary-value 
problem for buckling of doubly-curved shallow shells. The smallest value of 

the loading parameter p for which the equations are satisfied constitutes 
buckling of the shell. The equations are nondimens ional and are of order 
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one. Thus, the magnitude of the parameters multiplying the derivatives is a 
direct indication of the effect each term in the equations will have on the 
solution of a given problem. The parameters and Z^ are generalizations 

of the Batdorf Z parameter presented in references 11 through 13 for 
isotropic cylindrical panels. These generalizations of the Batdorf Z 
parameter rely heavily on the work presented by Stein in reference 8, and 
thus are referred to herein as the generalized Batdorf-Stein Z parameters. 

As the shell approaches a flat plate, the parameters Z^ and Z ^ approach zero 

and the equations uncouple. 

Parameters in Terms of Membrane Stiffnesses 

The parameters q^, fi , S^ 9 Z^ t and Z ^ have been given in terms of 

the coefficients of the inverted form of the membrane constitutive equation. 
Expressions for these parameters in terms of the usual membrane stiffness 
coefficients are obtained as follows. 

First, inverting the membrane stiffness matrix [A] associated with 
equations (35) yields 


11 ” ^ A 22 A 66' 

A 26 )/det(A) 

(59a) 

22 " (A 11 A 66‘ 

A* fi )/det(A) 

(59b) 

12 “ (A 16 A 26" 

A 12 A 66 )/det(A) 

(59c) 

66 “ (A 11 A 22’ 

A l 2 )/ det ( A ) 

(59d) 

l 16 “ (A 12 A 26' 

A 16 A 22 )/det(A) 

(59e) 

'26 “ (A 12 A 16' 

A ll A 26 )/det(A) 

(59f ) 
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where 


det(A) - (A 11 A 22 - A 12 )A 66 - A u A 26 - A 22 A lfi + 2A 12 A 16 A 26 ( 59 g) 

Substituting these expressions into the expressions for the nondimens ional 
parameters gives 

2 11/4 


m 


A 11 A 66 ' A 16 
2 

L A 22 A 66 ‘ A 26 j 


( 60 a) 


A 11 A 22~ A 12~ 2 A 12 A 66 + 2 A 16 A 26 


(A 11 A 66 ' A 16 )(A 22 A 66 ' A 26 ) 


1/2 


( 60 b) 


A 11 A 26 ~ A 12 A 16 


m 


2 3 


(A 11 A 66 ‘ A 16 ) (A 22 A 66 ' A 26 ) 


1/4 


( 60 c) 


and 


8 - 

m 


A 22 A 16 ~ A 12 A 26 


(A 11 A 66 ' A 16 )(A 22 A 66 ' A 26 ) ^ ^ 


( 60 d) 


Z, - 


2 

h 

( A H A 22 - a 12 )A 66 - A u A 26 - A 22 A 16 + 

L 6 A 26 

R i 

12 

(A 11 A 66 ' A 16 )(A 22 A 66 ' A 26 )D 11 D 22 

J 

1/2 

2 

L 2 

U. _J 

( A 11 A 22 - A 12 )A 66 ' A 11 A 26 ‘ A 22 A 16 + 2 A 12 A 16 A 2 el 

R 2 

12 

(A 11 A 66 ' A 16 )(A 22 A 66 ' A 26 )D 11 D 22 

1/2 


- 


_ 


1/2 


1/2 


( 60 e) 


( 60 f) 


For the case when . and k , r are zero-valued; 7 - 0 , 5 * 0 , and the 

I rj / r\ III TTl 


:>6 


m 


m 


remaining nonzero parameters simplify to 


m U 


2u 

L A 22. 


1/4 


( 60 g) 




A 11 A 22~ A 12 2A 12 A 66 


( 60 h) 


2A 66^ A 11 A 22 


26 



2 

, _h 

'1 R-. 


A 11 A 22' A 12 
127a 1] a 22 d 11 d 22J 


ll/2 


(601) 


"2 R„ 


A 11 A 22~ A 12 
127a 11 a 22 d 11 d 22J 


ll/2 


(60j) 


For the case of an isotropic material, the nonzero parameters are given by 


m L, 


M - 1 


z i - 5 ^ (1 - *' I)1/2 


(60k) 

(601) 

(60m) 


R 2 t 


(i - /) 1/2 


(60n) 


where t is the total plate thickness. 


Reduction to a Single Equation in Terms of W 

Generalized Donnell-Stein Equation . Donnell showed (see references 11 
through 13) that a single eighth-order differential equation could be 
obtained for isotropic cylindrical shells by eliminating the stress function 
appearing in the buckling equations. This task is performed by applying 
successive differentiation to obtain a single equation referred to by 
Batdorf as an escalated equation (see reference 14). Applying Donnell's 
approach to a doubly-curved shallow shell, equation (53) is operated on by 
the curvature operator to give 

D (D (SF) ) - yl 2 D 2 (6W) (61a) 

cm c 
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Using the fact the the differential operators in equation (61) are linear, 
their order of operation can be interchanged to give 

D m (D c (6F)) D c (fiW) (61b) 

Next, operating on equation (56) with the membrane stiffness operator yields 

W* W)) + ^ D m (D c (5F)) “ P D ra K g ( * W) <62) 

Substituting the right-hand-side of equation (61b) into equation (62) yields 
an eighth-order partial differential equation referred to herein as the 
generalized Donnell-Stein equation; i.e., 

D m (D b (SW)) + 12D*(SW) = pD m K g (SW) (63) 

Modified Batdorf -Stein Equation . The doubly-curved shallow shell 
counterpart to Batdorf' s equation is obtained by expressing equation (53) as 

6F - Jl2 D m 1 (D c (5W)) (64) 

where D ^ ( ) denotes the inverse differential operator (anti-derivatives) 

which symbolically represents integration. A detailed account of the use of 
inverse differential operators is given in references 11 through 14. Using 
the fact that the differential operator and inverse differential operator in 
equation (64) are linear, their order of operation can be interchanged to 
give 

SF - Jl2 D (D _1 (6W) ) (65) 

c m 

Substituting equation (65) into equation (56) yields the desired equation; 
i.e. , 

D, (6W) + 12D 2 [D' 1 (SW)] - pK («W) - 0 (66) 

b c m g 
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This integro-differential equation is a generalization of Batdorf's modified 
equilibrium equation presented in reference 12 and is referred to herein as 
the Modified Batdorf -Stein equation. 

Implied Membrane Boundary Conditions . As pointed out by Batdorf in 
references 11 through 13, the original Donnell and Modified Batdorf 
equations for isotropic cylindrical shells implicitly prescribe boundary 
conditions on the membrane displacements and stress resultants. This 
attribute of these equations is a direct consequence of the elimination of 
the stress resultant function. Thus, the Generalized Donnell-Stein and 
Modified Batdorf -Stein equations presented in this paper also possess this 
attribute. Results presented by Batdorf [11-13], Rehfield and Hallauer 
[15], and Sobel, et. al. [16] suggest that the effect of membrane boundary 
conditions on the predicted buckling load and mode shape of shallow 
isotropic cylindrical shells may, in several cases, be small. For these 
cases, buckling results obtained from the Generalized Donnell-Stein and 
Modified Batdorf - Stein equations give reasonable estimates of the collapse 
load. This information, however, is not generally known for shallow 
laminated composite shells. It does seem reasonable that, for the class of 
slightly curved shells which are not imperfection sensitive, the coupling of 
membrane and bending behavior is mild enough that the analytical predictions 
of buckling that are obtained from the Generalized Donnell-Stein and 
Modified Batdorf - Stein equations are reasonably accurate and thus useful in 
preliminary design of shallow shells. 

The actual membrane boundary conditions implied by solutions to the 
Generalized Donnell-Stein and Modified Batdorf -Stein equations are 
determined from the displacement form of the membrane equilibrium equations. 
These equations are obtained by substituting constitutive equations (35a) 
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through (35c) into membrane equilibrium equations (36a) and (36b). This 
step yields two coupled partial differential equations that relate the 
unknown membrane buckling displacements 8 u^ and 8 u^ to the known transverse 

buckling displacement 8 w. Using differentiation, the two coupled equations 
are converted to two independent uncoupled equations; i.e., one equation 
relating 8 u^ to 8w and another relating 8\x^ to £w. Performing these 

operations and nondimensionalizing the resulting equations yields the 
following equations for the special case of and being zero-valued. 


D m («V - - 

m 1 


a Z, + 
m 1 


12 


^ A ll A 22 


a 3 gw 

3 

dz. 


(2m + 


12 


7 a 1 1 A 2 2 


)Z, 


T z 2 

t 

m 


a 


2 

dz^dZp 


(67a) 


V'V 


12 1 

~Z===Z \ + 2 z 2 

/ A U A 22 “m 


3 3 5W 

3 

3z 2 


(2m + 


12 


y A u A 22 


> Z 9 ' Z 1 

z mi 


a 3 £w 

2 


(67b) 


where the nondimens ional membrane displacements 5U^ and 5 are given by 


6 U 


A 11 A 22~ A 12 


12^ A 11 a 22 D 11 D 22 


£u- 


(68a) 


and 


A 11 A 22' A 12 

5U - L 6 u 

12’/ A x 1 A 22 D 11 D 22 


(68b) 


These differential equations are generalizations of those presented by 
Donnell (see reference 11) for isotropic cylindrical shells. 
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A general method of determining the membrane boundary conditions 
implied by 6W is to expand both <5U^ and into full double Fourier series 

and then determine the Fourier coefficients using the orthogonality of the 
functions forming the double Fourier series. Once the Fourier series are 
determined, stress resultants are computed directly from the membrane 
constitutive equations. The boundary conditions are identified by examining 
the functional form of the Fourier series for the displacements and stress 
resultants on the boundary of the shell. 

Approximate Solution of the Buckling Equations 
Simple closed form solutions to equations (63) and (66) are not readily 
available since the geometric stiffness operator generally has variable 
coefficients and both the geometric stiffness and bending stiffness 
operators have odd-combination mixed partial derivatives (e.g., three 
derivatives with respect to and one with respect to z ^) . Approximate 

solutions, however, can be obtained by methods such as the Bubnov- Galerkin 
method (see reference 12 for example) . The usual way of solving equations 
(63) and (66) by the Bubnov- Galerkin method involves using series expansions 
for 8 F and 6W that satisfy all the boundary conditions of a given problem. 
These expressions are then substituted into equations (63) and (66) to 
obtain two residuals (one for each equation) since the series expansions 
generally do not satisfy the two differential equations. Then, the 
residuals are expanded in the same series (i.e. , same basis of the solution 
space) and the components of the residuals are forced to be zero-valued in 
an integrated sense. This process results in a generalized algebraic 
eigenvalue problem whose eigenvector includes the buckling displacement 
modal amplitudes and the stress function amplitudes. 
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Often in performing analysis for preliminary design, or parametric 
study, it is desirable to simplify the analysis as much as possible. One 
such simplification was introduced in shell buckling analysis by Batdorf 
(see references 12 and 13). The method associated with this simplification 
consists of a Bubnov- Galerkin solution to Batdorf' s counterpart of the 
Modified Batdorf -Stein equation. The results and discussion presented by 
Batdorf in references 12 and 13 suggest that use of the Modified Batdorf- 
Stein equation is preferred to the use of the Generalized Donnell-Stein 
equation when applying the Bubnov-Galerkin method. By eliminating the 
stress function using the inverse differential operator, the size of the 
generalized algebraic eigenvalue problem to be solved in the Bubnov-Galerkin 
method is significantly reduced, and thus the attractiveness of the analysis 
for parametric study is greatly increased. The application of this method 
to the Modified Batdorf -Stein equation is outlined in the subsequent 
section, 

Bubnov-Galerkin Formulation for the Modified Batdorf -Stein Equation 

To obtain an approximate solution to equation (66), the nondimens ional 
transverse buckling displacement is expressed as 

M N 

5 W " 5 W mn~ VpqWVV (69) 

p-1 q=l n 

where the basis functions $ (z. ) and * (z 0 ) are selected such that 6W 

p i q l MN 

satisfies all the boundary conditions of the given problem, and are at least 

a linearly independent set. Substituting equation (69) into equation (66) 

yields a residual since the assumed displacement series does not satisfy the 

integro-differential equation, equation (66). The residual is given by 
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( 70 ) 


M N 2 .1 

£ £ a {D, ($ * ) + 12D (D ($ * )) - pK ($¥)}- 0 


pll qll PI b ? 0 


cm p q 


g p q 


and Is dependent on the number of terms taken in the series. The evaluation 

of the terms involving D and K are straightforward, however, the 

^ g 


evaluation of the term involving the inverse operator is somewhat involved. 


The general procedure for determining D^($ ¥ ) is to first expand ^^(SW) 

in a general double Fourier series as if it is an arbitrary function of two 

variables (see reference 14). Next, the condition that D (D ^(6W)) — SW is 

mm 

enforced by operating on the Fourier series with (i.e., the left-hand- 

side of the equality) and then substituting equation (69) for 5W appearing 
on the right -hand- side of the equality. The Fourier coefficients are then 
determined using the orthogonality property of the functions forming the 
basis of the Fourier series. The resulting equations yield the following 
equation 


D 1 (6W) - 
m 


D 


00 

4 


M 

E 

p-1 


N 

e : 

q=l 


m 


($ * )a 

p q pq 


(71) 


The leading constant term is of no significance in the Bubnov-Galerkin 

2 -1 

solution since it vanishes once D (D (6W)) is computed. 

c m 

In the Bubnov-Galerkin method, the residual is assumed to be a function 
that can be expanded in the same basis functions as 6W ; i.e., it is assumed 
that, the residual can be expressed as 


R, 


MN 


M N 

E E r $ (z 1 )$ <z 0 ) 

p-l q-1 W P 1 q 2 


(72) 


The coefficients r represent the components of the residual in the 

PQ 


function space spanned by and . The coefficients r are determined by 
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multiplying equation (72) through by the set of functions (for m = 1, 

2, . M and n - 1, 2 , N) and integrating the resulting equation over 

the region on which the boundary -value problem is defined. This step gives 
i l 


RmN* 


m (z, )4 (z )dz dz 
mini Li 


M N 

S Sr 


p=l q=l 


pq 


F H 
pm qn 


(73) 


where 


and 


F 

pm 


i 

* (z, )$L (z n ) dz i 

m 1 p 1 1 

o 


l 


H 

qn 




(74a) 


(74b) 


An approximate solution to the boundary- eigenvalue problem is obtained by 
requiring that the components of the residual series expansion be zero- 
valued, Noting that the F and H coefficients of the right-hand-side of 

pm qn ° 

equation (73) constitute a nonsingular coefficient matrix of the linear 
system of equations defined by equation (72) , a sufficient condition for the 
components of the residual series to vanish is given by 

l l 

Wm< Z l>V Z l> dz l dz 2 - ° (75) 

J oJ 0 


for all combinations ofm-1, 2, . . . , M and n = 1, 2, . . . , N. More 

precisely, the following equations must be satisfied. 


M N 
S S K 


a 

mnpq pq 


M N 

S S K a 
p=l q=l mn P q Pq 


(76) 


p=l q=l 

where K constitutes the stiffness matrix given by 


mnpq 


K - K B + K C 
mnpq mnpq mnpq 


(77a) 
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The stiffness coefficients associated with bending are given by 


K 


B 

mnpq 


1 1 


D, ($ $ )$ dz-dz 0 
b p q m n 1 z 


(77b) 


and the stiffness coefficients associated with shell curvature are given by 


K 


mnpq 


12D Z (D' 1 <$ ))$ dz,dz 0 

cm pq m n 1 z 


(77c) 


The term K constitutes the geometric stiffness matrix and its 

mnpq 

coefficients are given by 


K 


mnpq 


K ($ $ )$ $ dz^dz 0 
g p q m n 1 z 


(77d) 


Once the specific form of $ and ^ are given, the inverse differential 


m 


n 


operator can be obtained, the stiffness and geometric stiffness matrices can 
be computed, and the eivenvalue problem can be solved. Buckling is defined 

by the smallest value of p that satisifes the generalized algebraic 
eigenvalue problem defined by equation (76). 


Results and Discussion 

Following the procedure presented in the previous section of this 
paper, the buckling behavior of a shell can be determined in terms of a^, 

7i , 8^ , a , u, 7 , 6 , Z, , and Z 0 . For a given family of laminates, such as 
b b m m m 1 z 

the f(+0) 1 laminates (n - 1, 2, ...) described in reference 1, changes in 

the fiber orientation and stacking sequence of a laminate generally results 
in changes in all the nondimens ional parameters presented herein. Plots 
showing the dependance of several of the parameters on fiber orientation and 
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stacking sequence are presented in reference 1. The point to be made is 
that the parameters are not independent with respect to laminate 
construction. However, insight into determining, and understanding, the key 
parameters affecting buckling behavior can be obtained by studying how the 
parameters vary with respect to laminate construction, and by studying the 
sensitivity of the buckling behavior with respect to varying each parameter 
in an independent manner. For example, an indication of the sensitivity of 
the buckling behavior with respect to variations in these parameters can be 
obtained from plots of a buckling coefficient as a function of each 
parameter. To demonstrate this philosophy, some typical results of a 
parametric study are presented in the next section of this paper for a 
representative example problem. 

Example Problem 

The example problem presented in this section is a shallow shell that 
is loaded on its edges by a uniform shearing traction r as shown in figure 
2. The shell is supported such that the transverse displacement and 
rotation along the edges are zero-valued (clamped with respect to bending 
behavior). The basis functions used in the approximate analysis presented 
herein that satisfy these boundary conditions are given by 

* cos (p- l)?rz^- cos(p+l)irz^ (78a) 

$ - cos (q- 1 )ttz 0 - cos (q+l)7rz 0 (78b) 

q 1 1 

for p - 1, 2, ..., M and q - 1, 2, N. Associated with these basis 

functions, and the use of the Modified Batdorf -Stein equation, are implied 
membrane boundary conditions. To greatly simplify the analysis, the example 
problem is defined as a shell that has = A^^ “ 0; i.e., its membrane 
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behavior is specially orthotropic. The corresponding implied boundary 


conditions are given by 6 - 0 on the edges « 0 and 1, and 3 

6Ni2 “ 0 on the other two edges. 

Applying the procedure for approximate solution of the Modified 
Batdorf -Stein equation described in the previous section of this paper 
yields explicit expressions for the stiffness and geometric stiffness 

g 

coefficients defined in equations ( 77 ). The stiffness coefficients K 

^ mnpq 

Q 

and geometric stiffness coefficients K are identical to those obtained 

a mnpq 

from a Bubnov-Galerkin analysis of a flat plate and are not presented 

herein. The contribution of the stiffness coefficients associated with 

shell curvature to equation ( 76 ) is expressed by 

, M N 

^7 S £ K a « Z(m- 1 , n- 1 ) a 0 0 -{Z(m-l,n-l) + Z(m+l,n-l))a 9 + 

4 mnpq pq m-Z,n-z m,n-z 

t r p=l q=l 


Z(m+1 , n- 1) a - 0 -{Z(m-l,n-l) + Z(m-l,n+l))a 9 + 

m+z , n-z m-z,n 

{Z(m-l,n-l) + Z(m-l,n+l) + Z(m+l,n-l) + Z(m+1 ,n+l) la^ ^ + 

- { Z(m+1 ,n- 1) + Z(m+l,n+l) }a m+2 n + Z(m-l,n+l) a m _ 2n+2 + 

- { Z(m- 1 , n+1) + Z(nH-l,n+l))a +2 + Z(ra+l,n+l) a m+2i n+2 (79) 


where 

12 [Z 2 (m-l) 2 + Z 1 (n- 1) 2 ] 2 

Z(m-l,n-l) - 2 i 2 2 4 2 4 

{ a (m-1) + 2/i(m-l) (n-1) + (n-1) /a }ir 
m ra 

for m — 1 , 2 , . . . , M and n- 1 , 2 , .... N. 


( 80 ) 
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Results obtained from the approximate analysis are presented in figures 
2 and 3. The buckling resistance of a shell is indicated in these figures 
by the nondimens ional shear buckling coefficient defined by 


K - 

s 


r L„ 


2 

7T (D 


d !„) 1/4 


11 22 


(81) 


where r is the applied shearing traction. For the results presented in the 
figures, the orthotropic parameters (a^, /9, and /i) and the anisotropic 

parameters ( 7 ^, 8^ t 7 , and 5^) that are not varied are set equal to one and 

zero, respectively. This baseline set of values corresponds to an isotropic 
shell with sides of equal length (L^ - L^) . 

Results showing shear buckling coefficient as a function of the 
Batdorf -Stein shell curvature parameters, and , are presented in figure 

2. Results are shown in this figure for flat plates and for shells with 
zero, negative, and positive Gaussian curvature with values of Z^ and Z^ 

ranging from 0 to 100. The results presented in figure 2 indicate that the 
shear buckling resistance of a shell is significantly influenced by shell 
curvature, especially for the larger values of Z^ and Z^ shown in the 

figure. The results also indicate that the shells with positive Gaussian 
curvature are the most buckling resistant. Moreover, the shells with 
positive Gaussian curvature are more buckling resistant than those with 
negative Gaussian curvature, which are more buckling resistant than those 
with zero Gaussian curvature. Flat plates exhibit the lowest buckling 
resistance . 

Results showing shear buckling coefficient as a function of the 
curvature parameters, Z^ and Z^ , and the bending anisotropy parameters 7 ^ 
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and $ b (see equations (45c) and (45d)) are presented in figure 3. Results 

are shown in this figure for flat plates and for shells with positive 
Gaussian curvature corresponding to ** - 100. Values of the 

anisotropic parameters range from 0 to 0.5. This range of values is 
considered to be representative of a large class of laminated plates [1]. 
Results are presented in figure 3 corresponding to both positive and 
negative directions of the applied shear traction. The distinction of 
loading direction results from the presence of the bending anisotropy. 

The results presented in figure 3 indicate the shear buckling 
resistance of a shell with positive Gaussian curvature is more sensitive to 
variations in the anisotropic parameters than a corresponding flat plate. 

The results show substantial reductions in buckling resistance with 
increasing values of the anisotropic parameters for shells loaded in 
positive shear, and similar increases in buckling resistance for shells 
loaded in negative shear. Similar results were obtained for a corresponding 
shell with negative Gaussian curvature that indicate the same trend, but not 
to as large an extent as exhibited by the shell with positive Gaussian 
curvature . 

The results presented in figures 2 and 3 show that varying parameters 
independently can give insight into the factors driving the structural 
response. For example, by independently varying the parameters associated 
with shell curvature it has been found that positive values of Gaussian 
curvature substantially improve the shear buckling resistance of a shell. 

In addition, it has been determined that shell curvature can significantly 
affect the importance of the bending anisotropy on the shear buckling 
resistance. Both of these observations clearly indicate the benefits of 
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using nondimensional parameters to formulate the analysis and to perform 
parametric studies . 


Concluding Remarks 

A method of deriving nondimensional equations and indentifying the 
fundamental parameters associated with bifurcation buckling of shallow 
shells subjected to combined loads has been presented. Analysis has been 
presented for symmetrically laminated doubly-curved shells that exhibit both 
membrane and bending anisotropy. The analysis includes equations for 
nonlinear deformations and buckling of thin elastic shallow shells, and the 
procedure and rationale used to obtain useful nondimensional forms of the 
transverse equilibrium and compatibility equations for buckling are 
discussed. Fundamental parameters of the problem have been identified that 
explicitly indicate, in a compact manner, how both membrane and bending 
orthotropy and anisotropy influence buckling behavior. Generalizations of 
the well-known Batdorf Z parameter for symmetrically laminated shells with 
full anisotropy have also been presented, as well as generalized forms of 
Donnell's and Batdorf' s equations for shell buckling. In addition, shell 
boundary conditions and approximate solution methods of the nondimensional 
boundary- value problem have been briefly discussed. 

Results obtained from a Bubnov-Galerkin solution of a representative 
example problem have also been presented. The results demonstrate the 
advantages of formulating the analysis in terms of nondimensional parameters 
and using them to perform parametric studies. The results specifically show 
that shells with positive Gaussian curvature are much more shear buckling 
resistant than corresponding flat plates and shells with negative and zero 
Gaussian curvature. In addition, the results show that the importance of 
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bending anisotropy on shear buckling resistance is affected by shell 
curvature . 
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Figure 2. Effect of shell curvature on shear buckling resistance (L^/L^-l). 
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Figure 3. Effects of anisotropic parameters and shell curvature on shear buckling resistance 
(Li/L 2 - 1) . The parameters and are nondimens ional bending stiffness 

parameters defined by equations (45c) and (45d) . 
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